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Abstract 

Systems coupling fluids and polymers are of great interest in many branches of sciences. One 
of the models to describe them is the FENE (Finite Extensible Nonlinear Elastic) dumbbell 
model. We prove global existence of weak solutions to the FENE dumbbell model of polymeric 
flows for a very general class of potentials. The main problem is the passage to the limit in 
a nonlinear term that has no obvious compactness properties. The proof uses many weak 
convergence techniques. In particular it is based on the control of the propagation of strong 
convergence of some well chosen quantity by studying a transport equation for its defect 
measure. 

1. INTRODUCTION 

Systems coupling fluids and polymers are of great interest in many branches of applied 
physics, chemistry and biology. They are of course used in many industrial and medical 
applications such as food processing, blood flows... Although a polymer molecule may be a 
very complicated object, there are simple theories to model it. One of these models is the 
FENE (Finite Extensible Nonlinear Elastic) dumbbell model. In this model, a polymer is 
idealized as an "elastic dumbbell" consisting of two "beads" joined by a spring which can be 
represented by a vector R (see Bird, Curtis, Amstrong and Hassager [3, [8], Doi and Edwards 
[18] for some physical introduction to the model and Ottinger [55] for a more mathematical 
treatment (in particular the stochastic point of view) of it and Owens and Phillips [57] for 
the computational aspect). In the FENE model ([I]), the polymer elongation R cannot exceed 
a limit Rq. This yields some nice mathematical problems near the boundary, namely when 
\R\ approaches Rq. At the level of the polymeric liquid, we get a system coupling the Navier- 
Stokes equation for the fluid velocity with a Fokker-Planck equation describing the evolution 
of the polymer density. This density depends on t, x and R. The coupling comes from an 
extra stress term in the fluid equation due to the microscopic effect of the polymers. This is 
the micro-macro interaction. There is also a drift term in the Fokker-Planck equation that 
depends on the spatial gradient of the velocity. This is a macro-micro term. The coupling 
satisfies the fact that the free-energy dissipates which is important from the physical point of 
view. Mathematically, this is also important to get uniform bounds and hence prove global 
existence of weak solutions. 

The system obtained attempts to describe the behavior of this complex mixture of polymers 
and fluid, and as such, it presents numerous challenges, simultaneously at the level of their 
derivation [15], the level of their numerical simulation [57\ 134] , the level of their physical 
properties (rheology) and that of their mathematical treatment (see references below) . In this 
paper we concentrate on the mathematical treatment and more precisely the global existence 
of weak solutions to the FENE dumbbell model ([I]). These solutions are the generalization 
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of the Leray weak solutions [43[ [42] of the incompressible Navier-Stokes system to the FENE 
model. 

An approximate closure of the linear Fokker-Planck equation reduces the description to 
a closed viscoelastic equation for the added stresses themselves. This leads to well-known 
non-Newtonian fluid models such as the Oldroyd B model or the FENE-P model (see for 
instance [19], [TO])- These models have been studied extensively. Guillope and Saut [26] [27] 
proved the existence of local strong solutions, Fernandez-Cara, Guillen and Ortega [22], [21] 
and [23] proved local well posedness in Sobolev spaces. In Chemin and Masmoudi [9] local 
and global well-posedness in critical Besov spaces was given. For global existence of weak 
solutions, we refer to Lions and Masmoudi [48J. We also mention Lin, Liu and Zhang [45J 
where a formulation based on the deformation tensor is used to study the Oldroyd-B model. 
Global existence for small data was also proved in |41[ [39] . 

At the micro-macro level, there are also several works. Indeed, from the mathematical point 
of view, the FENE model and some simplifications of it were studied by several authors. In 
particular Renardy [58] proved the local existence in Sobolev space where the potential U is 
given by U{R) = (1 — |i?| 2 ) 1_cr for some a > 1. W. E, Li and Zhang [20] proved local existence 
when R is taken in the whole space and under some growth condition on the potential. Also, 
Jourdain, Lelievre and Le Bris [33J proved local existence in the case b = 2k > 6 for a 
Couette flow by solving a stochastic differential equation (see also [31] for the use of entropy 
inequality methods to prove exponential convergence to equilibrium). Zhang and Zhang [61J 
proved local well-posedness for the FENE model when b > 76. Local well-posedness was also 
proved in [51] when b = 2k > (see also [36]). One of the main ingredients of |51] is the 
use of Hardy type inequalities to control the extra stress tensor by the H 1 norm in R which 
comes from the diffusion in R. In particular no regularity in R is necessary for the initial 
data. Moreover, Lin, Liu and Zhang [46J proved global existence near equilibrium under some 
restrictions on the potential (see also the related work [39]). Recently many other works dealt 
with different aspect of the system. In particular the problem in a thin film was considered in 
[111 , the problem of the long time behavior was considered in [601 130[ [T] . the problem of global 
existence in smooth spaces in 2D for some simplified models (when there is a bound on r in 
L°° ) was considered in [13[ 14?] [T4] [54"] , the problem of non-blow up criterion was considered in 
[40] . the problem of stationary solution was considered in [HI [TO], the study of the boundary 
condition at dB was considered in [28[ [50] . 

More related to this paper, the construction of global weak solutions for simplified models 
was considered in (3J HI [5j [62j [601 16] in the case the system is regularized by some diffusion 
in the space variable or by a microscopic cut-off. The case of the co-rotational model was 
considered in |49j . The co-rotational model preserves some of the compactness difficulties of 
the full model. It allows to get more integrability on the ip which makes the compactness 
analysis much simpler. 

We end this introduction by mentioning other micro-macro models. Indeed, a principle 
based on an energy dissipation balance was proposed in [12], where the regularity of nonlinear 
Fokker-Planck systems coupled with Stokes equations in 3D was also proved. In particular 
the Doi model (or Rigid model) was considered in |56j where the linear Fokker-Planck system 
is coupled with a stationary Stokes equations. The nonlinear Fokker-Planck equation driven 
by a time averaged Navier-Stokes system in 2D was studied in [13] (see also [H]). Recently, 
there were many review papers dealing with different mathematical aspects of these models 
[59] 1441 138] ■ In particular we refer to [38] for an exhaustive list of references dealing with the 
numerical point of view. 

1.1. The FENE model. A macro-molecule is idealized as an "elastic dumbbell" consisting 
of two "beads" joined by a spring which can be modeled by a vector R (see [8]). Before writing 
our main system ([1]), let us discuss the main physical assumptions that lead to it: 
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• The polymers are described by their density at each time t, position x and elongation 
R. This is a kinetic description of the polymers. 

• The inertia of the polymers is neglected and hence the sum of the forces applied on 
each polymer vanishes. We refer to [16] where inertia is taken into account. Moreover, 
the limit m goes to zero it studied where m is the mass of the beads. 

• The polymer solution is supposed to be dilute and hence the interaction between 
different polymers is neglected. This is why we get a linear Fokker-Planck equation. 
Let us also mention that there are models for polymer melts such as the reptation 
model (see for instance |55|). 

• The polymer is described by one vector R in B(0, Rq). Let us mention that there are 
models where each polymers is described by one vector R such that \R\ = 1 (the rigid 
case, see [14]) or by K vectors Ri, 1 < % < K (see [6]). Usually the difference between 
these models comes from the length of the polymers as well as their electric properties. 

• In the Fokker-Planck equation an upper-convected derivative is used. This is can be 
seen as the most physical one. Other used derivatives are the lower-convected and the 
co-rotational ones (see [3 [8]). The co-rotational one has the mathematical advantage 
that one has better a priori estimates (see |49j). 

• We neglect the diffusion in x in the Fokker-Planck equation. Indeed, this diffusion is 
much smaller than the diffusion in R. Actually, it makes the mathematical problem 
much simpler. 

Under these assumptions, the micro-macro approach consists in writing a coupled multi- 
scale system : 

dtu + (u ■ V)ii — vAu + Vp = divr, dfvu = 0, 



(1) 



d t ifj + u. = div R -VuRi(j + fiVip + VUtp 



Tij = f B (Ri g> VjU)ip(t, x, R)dR (VUip + f3Vip).n = on dB(0, R ). 

In (pTJ) , tf)(t,x,R) denotes the distribution function for the internal configuration and F(R) = 
VrU is the spring force which derives from a potential U and U(R) = — Mog(l — |-R| 2 /|i?o| 2 ) 
for some constant k > 0. Besides, /3 is related to the temperature of the system and v > is 
the viscosity of the fluid. In the sequel, we will take f3 = 1. 

Here, R is in a bounded ball B(0,Rq) of radius Rq which means that the extensibility of 
the polymers is finite and if!! where f2 is a bounded domain of R D where D > 2 or Q, = T D 
or n = R D . In the case fi has a boundary, we add the Dirichlet boundary condition u = 
on d£l. We have also to add a boundary condition to insure the conservation of ip, namely 
(—VuRip + VUip + f3Vip).n = on dB(0, Rq). The boundary condition on dB(0, Ro) insures 
the conservation of the polymer density and should be understood in the weak sense, namely 
for any function g(R) £ C 1 (i?), we have 



dR. 



(2) d t / gnpdR + u.V x / g^dR = - V R g 

Jb Jb Jb 

Notice in particular that it implies that ip = on 8B(0,Rq) and that if initially J ip(t = 
0,x,R)dR = 1, then for all t and x, we have J ip(t,x,R)dR = 1. We will see later an other 
way of understanding this singular boundary condition. 

When doing numerical simulation on the FENE model, it is usually better to think of the 
distribution function ip as the density of a random variable R which solves (see [55] ) 

(3) dR + u.VRdt = (VuR - V R U{R))dt + V2dW t 

where the stochastic process Wt is the standard Brownian motion in WL N and the additional 
stress tensor is given by the following expectation r = E(i?, <g> V jU). Of course, we may need 
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to add a boundary condition for ([3]) if R reaches the boundary of B. This is done by requiring 
that R stays in B (see |32j). Using this stochastic formulation has the advantage of replacing 
the second equation of (|2.ip which has 2D + 1 variables by ([3]). Of course one has to solve 
([3]) several times to get the expectation r which is the only information needed in the fluid 
equation. This strategy was used for instance by Keunings [35] (see also [23]) and by Ottinger 
[51] (see also [25]). 

In the sequel, we will only deal with the FENE model and we will take (3 = 1 and Rq = 1. 



2. Statement of the results 

This paper is devoted to the proof of global existence of free-energy weak solutions to 
the FENE model. The main difficulty of the construction is the passage to the limit in an 
approximate system in the nonlinear term Vu n ip n . Indeed, we only have a uniform bound on 
Vu n in L 2 ((0,T) x 0) and ip n in L°°((0,T) x fijL 1 ^)) for all T > and so assuming that 
u n and ip n converge weakly to u and ijj, it is not clear how to deduce that Vu n Tp n converges 
weakly to Vuip. 

Before mentioning our main result, let us recall that the construction of global weak solu- 
tions to simplified models was considered in [H [5j [60j S3 E2] . In particular in [4J a diffusion in 
the space variable in the if) equation is added. Mathematically this yields a bound on V^y^ in 
L 2 ((0, T) x n x B) and hence one can easily pass to the limit in the product Vu n ip n using the 
Lions- Aubin lemma. This extra diffusion term is physically justifiable but it is much smaller 
than the diffusion in the R variable and this is why we did not include it here. Recently, 
Barrett and Suli [B] extended their results to the case of bead-spring chain models where each 
polymer is described by K springs R 1 , 1 < i < K again with diffusion in the x variable. Also, 
in |49j . the co-rotational model was considered. It allowed us to get more a priori estimates 
on if) n , namely one can get that ip n is in all LP spaces. An argument based on propagation of 
compactness similar to the one used in |48j allowed us to conclude. 

Here, we consider the more physical model (pQ). The system (pQ) has to be complemented 
with an initial data u{t = 0) = no and ip(t = 0) = ipQ. 

Notice that (u = 0,^oo) where ipoo 

e -U(R) 

(4) Voc(i?) 



f B e-U( R ')dR> 

defines a stationary solution of (pQ). To state our result, we first impose some conditions on 
the initial data. We take uo(x) G L 2 (S7), div(uo) = and i/jq(x,R) > such that po(x) = 
f ipodR G L°°($7). Here po(x) is the initial density of polymers at the position x. We also 
assume the following entropy bound : — ^ G L\ogL{Vt x B ,dxpQ{x)tl) 00 dR)) 1 namely 
(5) 

]—\\Llo S L(nxB, P0 (xU oo dRdx)= / ( ^log ; ; h Ijpo^ipoodRdx < oo. 

POlpoo B ( ' W ' J JnxB Po4>oo P()V>oo PQ^oo 

Finally, we also assume the following L^ 2 L log 2 L bound, that we will call "log 2 " bound: 
(6) / ^ — —rj-dx < oo. 



n 



1 + 



Notice that interpolating ([6]) with the L°° bound on p$, we can deduce the LlogL bound ([5]). 

Theorem 2.1. Take a divergence free field uq(x) G L 2 (Q) and ipo(x, R) > such that po{x) = 
J ipodR G L°°(J7) and ([S]) and §6§ hold. Then, JT]) has a global weak solution (u,ip) such 

that u G L°°(1R + ;L 2 ) n ^(R^H 1 ), G L°°(1R + ; L\ogL(Q x B, dxp{x)il) 00 dR))) where 
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p(x) = J B ipdR and w G L 2 (R + ; L 2 (ft; H^tp^dR))) and ()32p holds with an inequality < 

instead of the equality and (|42p holds (with £1 replaced by any compact K of Q, in the whole 
space case). 

Remark 2.2. 1) Of course u and ip have also some time regularity in some negative Sobolev 
spaces in x and R. This allows to give a sense to the initial data (see [48 J for more details). 

2) By f G LlogL(Q x B,dxp(x)ipoadR) we mean that f f nxB {f log / — / + l)p(x)t/: 00 dR < 
oo. Notice that ([5]) does not really define a norm. One can of course define a norm using 
Orlicz spaces. However, we do not need to do it here. 

3) If the domain Q has finite measure (bounded domain or torus) then, the extra bound © 



1/2 



dx < oo. This extra bound on the initial data allows us to 



reduces to f n \j B fa log ^ 

prove the extra bound ()42|) on the solution. This is useful to get some sort of equi-integrability 
of the extra stress tensor. Of course this is a very mild extra assumption, but it would be nice 
to see if one can prove the same result without it. Moreover, due to the local character of the 
weak compactness proof, the assumption (|42p can be weakened by assuming the bound to hold 



1/2 



locally in space, namely f K f B fa log 2 ^ dx < oo for any compact set K of 0. 

4) For the simplicity of the presentation, the proof will be given in the case po(x) is constant 
equal 1 and 0, has finite measure. We will also indicate the necessary changes to be done in 
the general case. 

The paper is organized as follows. In the next section, we give some preliminaries where 
we prove some Hardy type inequalities. In section 01 we derive some a priori estimates for 
the full model (pQ). In particular we recall the free energy estimate as well as a new "log 2 " a 
priori estimate which is useful in controlling the transport of the defect measures. In section 
[5l we prove the main theorem 12.11 As is classical when proving global existence of weak 
solutions, the only none trivial part is the proof of the weak compactness of a sequence of 
global solutions satisfying the a priori estimates and we will only detail this part of the proof. 
In section [6l we present one way of approximating the system. In section [7] we present some 
concluding remarks and some open problems. 

3. Preliminaries 

3.1. Hardy type inequalities. The dissipation term in the free energy estimate (|32p mea- 
sures the distance between ifi and the equilibrium tp^. We would like to use that bound to 
control the extra stress tensor in L 2 . This will be done using the following Hardy |29j type 
inequality. 



Lemma 3.1. Ifk>l, then we have 
(7) 

For k > 0, we have 
(8) 

For —1 < /3 < k < 1, we have 



o J 






and more generally for all 7 > 
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< c 



Jo \ X 



1-f 

2 




Remark 3.2. Before giving the proof, let us mention that this lemma should be compared to 
the results of section 3.2 of [51]. In particular Proposition 3.1 was used to control the extra 
stress tensor. However, the main difference is that the results of section 3.2 of [51] are done 
in an L 2 frame work since we were dealing with strong solutions there, however the results of 
lemma \3.1\ are in an L 1 frame work since we only have a control on the free energy and its 
dissipation. 

Inequality ([7|) for k > 1 is just Hardy inequality. Notice that there is no requirement on 
the boundary data since k > 1. To prove it, we make the change of variable y = x l ~ k and 

h(y) = Hence, to prove flj[J), it is enough to prove that 

poo r 2 poo yp. 

(11) I yy<c\ x h'(yf + —dy 

where a = > 1. To prove (jll|) . we integrate by parts in 

f A hh' f A h 2 , h(A) 2 h(l) 2 

for each A>1. The left hand side is bounded by C(J A ^dy) 1 / 2 (J A h' (yfdy) 1 ! 2 . To bound, 
h(l) 2 by the right hand side of (fTTj) . we use that h{y) < Cy/y since h'(y) 2 dy < oo hence, 
^s- goes to zero when y goes to infinity. This yields that 

f°° / h 2 \' r°° h h 2 

(13) h 2 (l) = - r-)dy = - 2-^-ti - a-^—rdy 



i \vj ' Ji ir y a+1 



which is controlled by the right hand side of (llip using Cauchy-Schwarz and the fact that 
a > 1. Letting A go to infinity, we get the result. 

The proof of ([8]) when k > 1 follows by interpolation. 

In the case < k < 1, (|T|) only holds if we add a vanishing boundary condition at x = 0. 
However, we can still prove that ([8]) holds without any extra condition. Indeed, making the 

1 & fxxTVt qti If 1 "\ q n r\ r\ on nlin re hi (a i\ . / ^(^P 



change of variables y = x (when k < 1) and denoting h{y) = y *gr, we see that ([8]) is 
equivalent to 

(14) H 1 y a ~ l h 2 dy} 2 < C U 1 y 2 «h 2 dy^j h'(y) 2 + y 2 «h 2 

where a = jzx- To prove (|14p . we integrate by parts in the following integral : 



(15) J\«hh'dy = -^j\^h 2 + f ^±. 



and notice that the left hand side is bounded by (^f^ y 2a h 2 f^h^y) 2 ^ using Cauchy- 
Schwarz inequality. 
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Moreover, we have 

h(l) 2 = [\y 2a+1 h 2 )'dy = [ 1 y 2a+1 hh' + (2a + l)y 2a h 2 
Jo Jo 

(16) < c(j\'(y) 2 + y 2 «h 2 j\ 2a h 2 ^' 2 . 

Hence, (I14D follows. 

When k = 1, we make the change of variable y = — logx and hence ([8]) is equivalent to 

(pea \ 2 / pea \ / pea \ 

J e-yh 2 dy\ <C(J e~ 2 yh 2 dyj (j ti \yf + e^h 2 ) 



and the proof of <\17h can be done in a similar way as that of (114ft . 

To prove ([9]), we first notice that if — 1 < j3 < 0, then the inequality can be easily deduced 
from ([8]) by interpolation. When (3 > 0, ([9]) is equivalent (in the case k < 1) to 

l 2 l — l — 

(18) (J y a ?- l h 2 dy^ <c(J y^dy) 2 Qf //(y^ + y 2 ^ 2 ) 2 

where ap = ■y^f and q = Applying ([14"]) with a replaced by otp, we get 

(19) N 1 y^'^dy) < C (£ y 2 ^h 2 dy\ ' (Y h'{y) 2 + y 2a h 2 ^ ' 

Notice that we kept a in the last term instead of putting ag. Indeed, the last integral comes 
from the estimate of h 2 (l) and we can keep a = in (fT6j) . Now, we can apply (fT9l) replacing 
otp — 1 by 2ap and we get 

(20) (J 1 y 2a ?h 2 dy^ < C (j* y^+^h 2 dy^j ^ (Y h'(y) 2 + y 2a h 2 ^j ^ 

We can iterate this, replacing a — 1 by 2a^, 2(2o^ + 1), ... in (fT9|) till we get an index greater 
than 2a = Interpolating with the last inequality, yields ([9]). 

In the case k = 1, Q is equivalent to 

1-/3 1+/3 

(21) ( jT e-^-^/i 2 ^ < C7 QH e- 2 ^ 2 dy) " ( h'(y) 2 + e" 2 ^ 2 dy) * 

The proof of (|2ip is similar and is left to the reader. 

For the proof of (fTU|) . we use that it is equivalent (in the case k < 1) to 
(22) 



2 -i— p 

V^/^log^ 2 )^ <cf C y 2a h 2 \og^{h 2 )dy\ 2 ( h!{y) 2 + y 2 ^ 



1-/3 . , . 1W 

27 



Again, one can prove (|22p in the case f3 = by an integration by parts similar to the one used 
in (|14p . The case where — 1 < (3 < can be deduced by interpolation from the case j3 = 
and the case < [3 < k can be deduced by a bootstrap argument as the one used in the proof 
of ©. 

3.2. Control of the stress tensor. We recall that Voo (R) = j t-u^) dR , = (l-|#| 2 ) fe//3 / / B (l- 
|^/|2-jfc//3 ^/ an j s i nce /5 = l ; ip^R) behaves like (1 — \R\) k when goes to 1. In particular 
we will apply lemma [37X1 with x = 1 — \R\. 

Using the inequality ([8]) in the radial variable with x = 1 — \R\, we get 
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Corollary 3.3. There exists a constant C such that we have the following bound 

I 2 



(23) 



|t(V)| 2 < ( / 1>dR) 



B 



B 



V 



R\ 



ipoodR 



This Corollary can be seen as the L 1 version of Proposition 3.1 of [51]. It will allow us to 
control the extra stress tensor by the free energy dissipation. 



3.3. Weighted Sobolev inequality. In subsection (J5JJ), we have to prove the equi-integrability 



of N£ ■ This will require the control of some higher LP space of y ^j—. We have the following 
proposition 

Proposition 3.4. There exists p > 2 and a constant C such that we have the following bound 



(24) 



B 



1 V 
V'oo 



X/p 



< 



1/2 



ipoo+i^dR 



For the proof we first notice that the only difficulty comes from the weight and hence we 
can restrict to the region where \R\ > ^. We also use some spherical coordinates, namely 
R = (1 — x)uj where to G S D_1 and < x < |. The square of the right hand side of (|24l) can 
be written as the sum of a radial part and an angular part : 



(25) 



(26) 



— I 2 x k dx duj. 




dto x dx. 



We recall the following ID weighted LP — L q Hardy inequality (one can also call it weighted 
Sobolev inequality) 



(27) 



1/2 
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F{x)\ q x k dx 



< C 



-1/2 



1/2 



\F'{x)\ 2 x k dx) . 

10 J \J0 J 

This inequality can be easily deduced from Theorem 6 of [37], taking u(x) = v(x) = x k for 

i(k+r 

k-l 

holds for any F, with F(^) = if 



any q < oo if k < 1 and for q < 2 [ k+ ^ if k > 1. Indeed, Theorem 6 of |37j stated that ([27 



0<r<-§ JO 



sup ( / x k dx) 1/q { / {x k )- l dx) 1/2 < 



oo. 



Hence, we get a control of J ^ in the space L 2 (§ D 1 ; ^^((0, |), x k dx)) using the radial part 
of the norm (|25p . 

On the other hand we can use the classical Sobolev inequality in D — 1 dimension to control 

in the space L 2 X ((0, ^);L s (8 £) - 1 ),x k dx) where s = if D > 3, s < oo if D = 3 

and s < oo if D = 2. Interpolating between the two spaces L^Ll and L 2 L^, we deduce the 
existence of some p > 2 such that ()24|) holds. 
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3.4. Young measures and Chacon limit. We recall here two important weak convergence 
objects used in this paper, namely the Young measure and the Chacon's biting lemma. Ac- 
tually, these two notions are very related as was observed in Ball and Murat [2j. 

Proposition 3.5. (Young measures) If f n is a sequence of functions bounded in L 1 (C/;lR m ') 
where U is an open set ofM. N , then there exists a family (vx)xeU of probability measures on 
M. m (the Young measures), depending measurably on x and a subsequence also denoted f n such 
that if g : M. m — >• R is continuous, if A <ZU is measurable and 

g(f n ) -» z(x) weakly in L\A;R), 

then g{.) G ^(W 71 ; v x ) for a.e. x G A and 

z[x) = \ g{\)dv x {\) a.e. x G A. 

In the case where f n is bounded in LP{U\ M m ) for some p > 1 (or when f n is equi-integrable) , 
we can always take A = U and we have (extracting a subsequence) 

g(f n ) ~* I gWdv x {\). 

Proposition 3.6. ( Chacon limit) If f n is a sequence of functions bounded in L l (U; M m ) where 
U is an open set ofW N , then there exists a function f G L 1 (?7;lR m ) ; a subsequence f n and a 
non-increasing sequence of measurable sets of U with linifc_Kx) C^iEk) = (where Cn is 
the Lebesgue measure on M. N ) such that for all k G N, f n — 1 / weakly in L 1 (C7 — E^W 11 ) as 
n goes to infinity, f is called the Chacon limit of f n . 

It is easy to see that if f n is equi-integrable then the Chacon limit of f n is equal to the 
weak limit of f n in the sense of distribution. 

If we consider continuous functions gt ■ M m —> W m , k G N satisfying the conditions : 

(a) <7fc(A) — > A when k — > oo, for each A G M m , 

(b) \g k (X)\ < C(l + |A|), for all k G N and A G M m , 

(c) limiAi^^ |A| -1 |5fc(A)| = for each k, 

then, under the hypotheses of Proposition 13.51 for each fixed k, the sequence of functions 
9k{f n ) is equi-integrable and hence (extracting a subsequence) converges weakly in L 1 (C/;M m ), 
to some fk- Applying a diagonal process, as k goes to infinity, the sequence fk converges 
strongly to some / in L 1 (C7;lR m ). The limit / is the Chacon's limit of the subsequence f n 
and it is given by 

f(x) = / \du x (\) a.e. x G U. 

This gives an other possible definition of Chacon's limit which is equivalent to the one given 
in Proposition 13.61 For the proof of these results we refer to [2]. 

4. A PRIORI ESTIMATES 

4.1. Free energy. The second equation of (pQ) can be written as 



(28) d t tp + u.Vip = divR 



Vu ■ Rip 



+ divR 



Woo 



We define p(t, x) = J B ipdR. Integrating (f28|) in R, we get the transport of p, namely dtp + 
u.Vp = 0. 

Multiplying ([28]) by log and integrating in R and x, we get 



(29) d t I I ^log(-^-)-^ + ^oo= I [ Vu-RV r U<iIj-A I [ ^ 

JflJB PWoo JnJB JVtJB 




dR 
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where we have used that V^oo = —ijjooVU. 

The first equation of (pQ) yields the classical energy estimate for the Navier-Stokes equation 



(30) 



d. 



= - Vu : t - v I \Vu\ 2 . 
Adding (j29|) and (f30j) yields the following decay of the free-energy 



(31) dtff [Vdog(^) - V + W>oo] 
Jn Jb PWoc 



-f I |Vu| 
In 



4 / / Voo 

In Jb 



V, 



Integrating in time, we get the following uniform bound for all i > 



Woo 



(32) / / [^log(-^)-^+^ 00 ] + ^ (*)+/"* v [ \Vu\ 2 +4 [ f Vc 
Jn Jb PVoo * Jo Jn Jn Jb 




C . 



To simplify the notations in the rest of this section, we will assume that po(x) = 1. The 
proof in the general case is identical and we will indicate the changes to be made at the end. 
The general idea is the following: When proving a priori estimates, one has just to replace 
Woo by p(t, x)if)oo and take advantage of the fact that p is just transported by the flow. When 
proving weak compactness, one can use that p n converges strongly to p in all L p ((0,T) x Q) 
spaces and use p n (t, x)V>oo- Due to the local character of the proof of weak compactness, a 
simpler way is to just use ■i/'oo and so the calculations given in section [5] hold even when po is 
not constant. 

4.2. log 2 estimate. The free energy only gives an LlogL^tp^dR) bound on For some 
technical reasons, we will need to control a slightly higher growth of tp in the R variable. 

We introduce %j) = ip + aV>oo for some a > 1 . This is done to insure that log ^— does not 
take negative values. It will also add a new term in the equation which will not present any 
extra difficulties. Hence, solves 



(33) 



dtW + u.S/ip = divR 



Vu-Rtp 



+ divR 



W ' 

Woo^R- 

Woo- 



aVu ■ RijjooV rU. 



We first derive this extra bound in the case the domain Vt is bounded and then discuss the 
modification of the argument in the whole space case. 

4.2.1. Case of a bounded domain. Multiplying (133j) by log 2 ^— and integrating by parts in R, 



we get 



W 



w 



(d t + u.V x ) / #og^(^-)-21og(-H + 2] = -2akV iU 



B 



Wo 



B 



RiRj 
1 - \R\ 



■Woo log 



(34) 



Vu-Ry2 log(— )- r V R — 

Li Woo V 



Wo 



B 



The second term on the right hand side of ()34p can be rewritten 



Woo 



Woo 



log( 



Woo 



Vu-R Voo V R [^\og{^)-^-) = 2 [ Vu-RV R U^ (log(-^)-l 

B \ Woo Woo Woo J B \ Woo 



w , , w 

Woo Woo W 

Taking the square root of ()34p . we get 
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(<;, + „.V,)| / #og 2 (-^)-21og(-^) + 2] 

IB Woo Woo 



1/2 



-akViUj J B jz^Woo lo 



/B#og 2 (£)-21og(i) + 2] 



(35) 



+- 



J B V,-^21og(£)^V R £ 

1/2 



I B Woo 



1/2 



log(#:) 



B 



#°g 2 (^)-2iog(^) + 2] 



/ B #og 2 (^)-21og(£) + 2] 



1/2 



h + h + h- 



Let us introduce the notation 



(36) 



N, 



#og 2 (-^)-21og(-^) + 2] 

B ^oo Woo 



1/2 



To bound Ii we use that, ib^ log 2 -±- < Ctb. Hence, the numerator of I± is bounded by 

C\Vu\ J j^pdR which is clearly in L 1 ^^) x 12 x B). Indeed, by using ([8]) and Corollary 
(|3.3p . we see that 



(37) 



Woo 



1/2 



To bound the second term on the right hand side of (|35p . we use that the numerator can 
be bounded by 



V7 DM , W ^Woor-j W 

Vu ■ Rip log(— )^Vi?— 

B Woo 1p Woo 



W) £ r\Y„\ j / Vool iog(-^-; 

'B Woo 



< 



Vb1 



V'oo 



2\ 1/2 



B Woo Woo 



1/2 



(39) < C|Vu| 2 ( / $log(-^)| I + ( / VUlogA 

\Jb Woo \Jb Woo 



^ V Woo 



(40) < C'|Vn| 2 (l + a) 1 /2^^ lo g2 ( _|_ ) ^ + ( J^°o| log(^-) 
Dividing (j38j) by A^2, we deduce that 



(41) 



/ 2 < C\Vu\ 2 - -h. 



Integrating (|35p in time and space and using the fact that I3 has a sign, we deduce the 
following a priori bound 
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(42) 



J s #og 2 (£)-21og(£) + 2] 



1/2 



< c T 



1/2 



for < t < T, if the initial condition satisfies f n [ f B tp [log 2 {^) - 21og(^) + 2] J < C . 

Hence, we see that (I35D can be written as 

(43) (d t + u.V)N 2 = F 2 

where F 2 is in L x ((0,r) x O). 

It turns out that passing to the limit in the bound (|42p is not clear. Actually, one can find 
sequences of functions ip n such that (|42|) holds and the weak limit does not satisfy (|42|) . This 
is the reason, we prefer to write the second bound as 



(44) 



B 



1/2 

g z logig^oodR ) dx + 



r 



+ 2] 



1/2 



where g is given by g = \J ^ log 1/2 (^). 



4.2.2. Case o/ unbounded domain. In the case = ~R D , we first take ci and C2 the two 
constants such that the function <j){x) = x[log 2 x — 2 log a; + c±] + c 2 satisfies the fact that 
cj)(l + a) = 4>'(1 + a) = 0. This is achieved by taking c\ = 2 — log 2 (l + a) and c 2 = 
2(1 + et)[log(l + a) — 1]. Notice also that the function (j>(x) is nonnegative for x > a since a is 
taken big enough. It is clear that the extra bound ([6]) implies that 



(45) 



/ 

Ju 



B 



i + 



1/2 



dx < a 







and hence, we can perform the same calculations as (i34|) and (f35|) with J B ■0[log 2 (^) — 

21og(^) + 2]dR replaced by j B 4>{-^)dR and with the function s — > y/s used to go from (j34)) 

which behaves like (f>i(s) = min(y / s, s). The rest of the proof is 



to (|35|) replaced by s 
identical. 



l+v^ 



4.2.3. Case p is not constant. In the case p is not constant and we are in a bounded domain, we 



have to modify (|34p slightly and multiply by log 2 



p4>° 



In the case we are also in an unbounded 



domain, we have to replace J B V'flog 2 



21 °g(#-) + 2]d J Rby/ J 



B 



■ (l+a)j> 



)dR. The extra 

factor is used to insure that when p is at microscopic equilibrium, namely ?/> = (p + a)V>oo, 
the integrand reduces to 0(1 + a). The rest of the proof is identical and yields at the end the 
following bound instead of (|42p 



(46) 




(<) + 



lrW . (1+a) ^ 
lQ g I (p+a) Voo 



1 + 



(I 



(1+a)^ 
B ^ (p+a)i/);x; 



1/2 



< Ct- 
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One can then deduce from (|46p that ()42|) and (|44p hold with the integration set O replaced 
by any compact K of MP . 



5. Weak compactness 

As it is classical when proving global existence of weak solutions, it is enough to prove 
the weak compactness of a sequence of weak solutions satisfying the a priori estimates of 
the previous section. In the next section, we present one way of approximating the system. 
We consider (u n ,ip n ) a sequence of weak solutions to (pQ) satisfying, uniformly in n, the free 
energy bound ([32]) and the log 2 bound (j4~2j) with an initial data (wqjV'o) sucn that (uQ,tpo) 
converge strongly to (uo,ipo) in I? (0) x Lj oc (£l; L 1 (£?)) and ipo log — ipo + p^tpoo converges 

strongly to ipo log ipo + poipoo in L l {VL x B). We also assume that (u n ,ij) n ) has some 

extra regularity with bounds that depend on n such that we can perform all the following 
calculations. 

We extract a subsequence such that u n converges weakly to u in L p ((0, T); L 2 (Q)) D 
L 2 ((0,T);H^(n)) and ip n converges weakly to ip in L p ((0, T); L] oc (Q x B)) for each p < oo. 
We would like to prove that (u,ip) is still a solution of ((!]). The main difficulty is to pass to 
the limit in the nonlinear term \7u n Rip n appearing in the second equation of ([T]). 

We introduce g n = J log 1 ^ 2 (^) and f n = \J ^ where tp n = ip n + a-^oo and a > 1 is 
any constant. We also assume, extracting a subsequence if necessary, that g n and /" converge 
weakly to some g and / in LP((0,T); Lf oc (tt x 5, dxipoodR)) for each p < oo. To prove that 

(u, ip) is a solution of (P), it will be enough to prove that (g n ) 2 = ^ log(^) converges weakly 
to g 2 = log(^-), namely that g n converges strongly to g in L 2 ((0, T); L 2 (QxB, dxipoodR)). 

First, it is clear that u, if) and g satisfy the same a priori estimates that the sequence u n , ip n 
and g" satisfy since all those functionals have good convexity properties. In particular it is 
clear that u, ip satisfy (|32p . We just point out that to pass to the limit in the last term on the 
left hand side of (j32|) . we can use the fact that the function (p2(x,y) = ^- is convex. To pass 
to the limit in (|4"4"|) . we also use the fact that 4>2(x,y) is convex. Hence, we deduce that 

(47) sup / ( ( [ g 2 \og{g 2 )^dR) 1/2 + Ni) dx(t) + f [ [ < C T 

o<t<TJn\ y JB J J Jo JnJB No- 



where N£ is the weak limit of ( f B ^ n [log 2 (|^) - 2 log (|^) + 2] 

Dividing (j33j) by tpoo we get 
(48) 



1/2 



ip r - 



oo 



ijj n ip n ip n 



d t -—+u n .V^— = div R -Vu n -R-— +VU.Vu n R— +A R ^ W.V^-f aVu n -R.V R U. 



Woo tpoo Wo 



From (|48p . we deduce that for any smooth function O from (0, oo) to M, we have 

J.n J.n Jjn J.n J.n 

dM^-) + « n .ve(£-) = -vu»fl • v*e(£-) + VrU.VuR^-q'(^) 

YOO YOO YOO YOO YOO 

(49) + A0 A - V R U ■ VrQ(^) - Q'\^)\V R fL\ 2 

Woo tpoo tpoo tpoo 

2akV u n RlRj Q'(^-) 
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We take 0(i) = t 1 / 2 log 1 / 2 (i) and recall that 5™ = 0(J^)- We introduce the following 
defect measures 7^ , 7^ • and such that 



VnV 1 -> Vug + 7, Vn™^e'(^) -> Vu^Q'(^) + 7' 

(50) Voo V'oo ^00 V^oo 

Vu n tjj n -> Vmp + 

where 7,7' G L 2 ((0,T) x O x fi) and /3 G £ 2 ((0,T) x tt;L l (B)) are matrix valued. 
On one hand, passing to the limit in (|49p with @(t) = t 1 / 2 log 1 / 2 (t), we get 



d t g + u.Vg = div R 



ViUjRjg - JijRj 



+ V R U.VuR—Q'{—) + V r UR : 7 

Woo Woo 



(51) 



1 

Woe 



Wooing 



+ 



|Vfl/n| 2 (logV 2 + log- 3 /2 )( ^) 



ak&(-j—)Viv2 J 



'Wo 



l-\R\< 



where, here and below, F n denotes the weak limit of F n and where we have used that 



(52) 



e / ( s ) = ^- 1 / 2 (iog 1 / 2 ( s )+iog- 1 / 2 (s )) 
e» = -i s - 3 / 2 (io g 1 /2( s ) + io g - 3 / 2 ( s )). 



Multiplying by g, we get 



d t g 2 + u.Vg 2 = div R -VuRg 2 + VuR ■ V R U[^-&{^-) )2g 

J V^OO Won ' 



divRfrijRjfig + VrUR : 7 '2 5 



(53) 



+ -— divR 

Woo 



Woo^R9 Z 



2|V B5 | S 



+ 



IV^PCW^ + log- 3 / 2 )^!)! 2 



/' 



2g -ak Q'(J—)V iU n- 

Won I 



n ^RiRj_ 
\R\- 



Multiplying (|53j) by V'oo and integrating in yields 



(d t + u.V) / VooS 



Woo (g 2 -2g^—Q'(^—) 

V Woo Woo 



(54) 



div R (7 i ji? J )25t-i/ , oo + VrUR : 7 2^ c 



+ I Woo — — 2g - 2y OQ \V R g\ 2 



WooakQ'(p)V lUl ^- 2 

B Woo J l-\R\ 2 
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where we recall that Tij(w) = 2k j B w ^zrm2 dR. Here, there is a small problem of definition: 
The terms on the second line of the right hand side are not well defined in the sense of 
distribution and we need some further analysis to make sense of them. Also, the transport 
term is not well defined even if we write it as div(u J B Wood 2 )- Actually, as we will see later, 
we will not use (|54p but a renormalized form of it. Indeed, we will construct in the next 
subsection a renormalizing factor N that satisfies (dt + u.V)-^ = and we will make sense of 
(|54p after dividing each term by iV 4 . 

On the other hand, passing to the limit in the equation satisfied by Wm we get 



(55) dti> + u.Vtp = divR 



Vu- Rip- fcjRj 



+ divR 



Woo^R 



Woo 



T-7 R{Rj 

2ak\/u : ipo 



1-\R\< 



Besides, ?/> n log(£) satisfies 



(d t + u n .v) 



(56) 



riog(^-: 

Woo 



4 / Woo 
' B 



V K 1 



' w r ' 
Wo 



Vu n : r(w n ) 



2ak I Va"-^%c, log! - 
B 1-\R\ ^o 



We would like to pass to the limit weakly in (|56|) and deduce that 



(d t + U .V) I ^log(^) 

IB Woo 



Vu : tW) + / ft, RlRj 



B 



l-\R 2 \ 



(57) 



Wc 



B 



' 1p r 

Wo 



2ak j 1 ( ,g(_)V^ i ^ Wo 



W T 
'Woo 



However, we can not use (|50p to pass to the limit in Vu n : r(^ n ) = J B Vu n ^f^w n and 
deduce that 



m Vu" :r(r" ] -Va:r(r ] + IJiJT^^\ 



since Vu n V*" i s om y bounded in L l (dtdxdR). Besides, we can not pass to the limit in 
the transport term even if we write it in divergence form. 

To overcome these difficulties, we will divide (|56[) by 1 + 8N2 where N% solves (|43p before 
passing to the limit. Then, we will send 5 to zero. 

To be able to deal with the limit 5 to zero, we need to renormalize (j56[) too. We denote 

A? = f B r log(£), iV 2 " = (j B #>g 2 (£) - 2 log(£) + 2]) and 6 K (s) = We first 

multiply ([56]) by 6' K (N™) and get an equation for 6 K (N™). Dividing the resulting equation by 
I+5N2, using (|43p and taking the weak limit when n goes to infinity (extracting a subsequence 
if necessary), we get for k, 5 > 
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{d t + U.V) 



(59) 



i + sm 



--Vu n 



T(tjj r 



(1 + 6N$)(1 + kN?) 2 (1 + 5N™)(1 + kN?) 2 Jb 



2a/c 



(! + <5iV 2 n )(l + KiV 1 n ) 2 



B V'oo 



(1 + <5iV 2 n ) 2 



?«(JVT). 



V'oo 



Now, we can send 8 to zero. Notice that due to the fact that 9 K (N™) is bounded and that 
F n is bounded in L , we deduce that the last term goes to zero when 5 goes to zero. Then, 
we send n to zero and recover at the limit 





= S,K f 

(d t + u.V)9 =Vu n : r(ip n ) - / Voo 

Jb 




1 Woo 


(60) 


"/. W, i-W 


»10g 


In 

7^ 

Woo 



t6,k 



-8,k 



where 9 



hm K -> hm^o i + gjv~ 



lim K _>.o 9k(Ni) is the Chacon limit of iV™ and 



-Fn* 5 K = h m h m 



«4oal>6 (1 + SN^)(1 + fdV™) 2 

for any sequence F n bounded in L 1 . We will prove in the next subsection that iV™ is equiinte- 
grable and hence that 9 is also the weak limit of . We also notice that if F n is equi-integrable 

then F n ,K = F n . 

To be really precise the term uXJ9 = div(#u) on the left hand side of (j60|) is not well defined 
since 9 is only in L^L\ and u is in L 2 H^. To give sense to (|60|) we need to use a renormalizing 
factor. Actually, we made the same remark for the transport term in the equation (|54|) . Recall 
that at the end, we would like to take the difference between (|6Up and (|54|) after dividing it 
by the factor iV 4 that we are going to define now. The fact of dividing by N will insure that 

is bounded and that all the terms will make sense. So the point is to divide (1591) by iV 4 
and then send 5 and k to zero. 

5.1. The renormalizing factor N. We recall that from, (j4~3|) . (3m(N%) solves 

(61) (d t + n".V)/3 M (iV 2 n ) = P' M (N2)F n . 
where f3u(s) = 9i/ M (s) = Passing to the limit in (jBTj) . we get 

(62) (d t + u.V)(3 M (N%) = P M {N2)F". 

This equation does not seem to be very useful since the right hand side is a measure. To 
overcome this problem, we first introduce the unique a.e flow X n in the sense of DiPerna and 
Lions |T7] of u n , solution of 



(63) 



d t X n (t, x) = u n (t, X(t, x)) X n (t = 0,x) = x. 



We also denote by X the a.e flow of u. 
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Let Q n be the solution of (|43p with F n replaced by \F n \ and taking the same initial data 
as N% at t = 0. Hence, 

(64) **^= A ,«T"(<,r(«,))l 

where the equation holds in the sense of distribution. 
Passing to the limit weakly in (|64|) . we get 



d[/3 M (Q n ){t,X"{t,x))} 



(65) ^"^ ^ = P' M m\F n (t,X n (t,x))\ 

From the stability of the notion of a.e flow with respect to the weak limit of u n to u, 
we know that X n (t,x) converges to X(t,x) in Lj and also that (X n (t)~ 1 )(x) converges to 
(X(t)~ 1 )(x) in Lj oc . This allows us to get the following equality of the weak limits 



(66) [p M (Q n )(t,X n (t, x))] = [P M (Q n )(t,X(t,x))}. 
Now, sending M to infinity in (|65p . we deduce that 

(67) d[Q(t,X(t,x))] = F 

dt 

where Q = limM-j-oc [PuiQ 11 )] is the Chacon limit of Q n and F = limM->oo (3' M (Q n )\F n (t, X n (t, x))\ . 
It is easy to see that Q G L°°(0, T; L 1 (0)) and that F G M({0,T) x fy. Integrating in t, we 
deduce that a.e in x G fi, we have 

(68) Q(t,X(t,x)) = Q{0,x) + f F(s,X(s,x))ds 



Jo 

for a.e t G (0,T). Due to the fact that F is nonnegative, we deduce that Q(t,X(t, x)) is 
increasing in time. We define the normalizing factor N by the following 

(69) N(t, X(t, x)) = Q(T , X(T ,x)) = Q(0, x) + f ° F(s, x)) ds 

Jo 

for t G (0, To) where To < T is a fixed time. In the sequel, we will denote T = To and 
will not make the distinction between these two times. Notice that ./V is constant along the 
characteristics of u, that N is in L°°(0, T; L 1 ^)) and that N(t,X(t,x)) is in L 1 ($7; L°°(0, T)). 
Moreover it is bounded from below by 1. Hence, it solves 

(d t + u.V)^ = 0. 



Also, the following two inequalities hold 



(70) Pm(N%) < MQ n ) <Q<N 

and hence the weak limit of N£ which is equal to the Chacon limit of N£ is bounded by N. 
The fact that the weak limit of N£ is equal to its Chacon limit comes from the fact that 
the sequence N£ is equiintegrable. This is a simple consequence of the dissipation of the free 

energy and the weighted Sobolev inequality (|24l) . Indeed, from (|24l) . we deduce that J ^ is 

bounded in L 2 ((0,T) x Q- LP (ijjoodR)) on the other hand from the conservation of mass, we 

know that \ l ^— is bounded in L°°((0,T) x Q; ^(ip^dR)). Interpolating between these two 

bounds, we easily deduce that */ ^— is bounded in T r ((0,T) x x B ^dtdxipoodR) for some 
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r > 2 and hence iV^ is equiintegrable. We also get that iVJ 1 is equiintegrable and hence 9 
which is the Chacon limit of iV™ is equal to the weak limit of iV". 



-5,K 



5.2. The term V« n : r(^ n ) . In this subsection, we will prove that Vu n : r(^ r 
T Ib fiv i~\R 2 \ ' This wm follow from the following two lemmas 

Lemma 5.1. 

j,k RiRj 



-<5, ft 



Vu : 



(71) 
where z 5)K 



Vn™ : r(^ n ) 
(1 + <5iV 2 n )(l + kN?) 2 " 7 B ~ l-|i2 2 



(l+<5AfJ)(l+KA r ") 2 



Lemma 5.2. z <5 ' K converges strongly to Vu n, ijj n = Vutp + j3 in L 1 ((0, T) x $7 x B; dtdx x-\r\ . 
when 5 goes to zero and then k goes to zero. 

Denoting T n > 5 > K = ^^1^1^^ , we get that 
Corollary 5.3. 

fi^ K r RR- 

(72) Vn™ : r(^ n ) ' = lim lim Vu n : r n ' 5 > K = Vu : r(ip) + / l —±-dR. 

k->0i5->0 Jb 1 — \R z \ 

Proof of LemmaUm The proof of <JZTJ> follows from the fact that z n > s < K Y^il 



is equi-integrable in L 1 ((0, T)x£lx B;dt dx ^rnP ) for (5, k fixed. Indeed, consider the real val- 
ued function $(x) = xlog(l+x)+l. It is enough to prove that <i)(z n ' <5A ) = $ f v "'''"' 

is bounded in X = L 1 ((0, T) x f] x B;dtdx^^). To simplify notation, we denote iV r 
(1 + <5iV^)(l + fdV™) 2 . Hence, it is enough to bound 



^(I+M^XI+kAT^ 

ider the real val- 
\ju n j> n 

^ oc (l+5A r ?)(H-KA r r) 2 



(73) 



+ — log 



in X (see definition above). 

To bound the first term appearing in (f73|) we use the Hardy type inequality ([8]) to get that 



(74) 
(75) 
(76) 



N n J B ^ S Uoo ) 1 - \R\ 



dR 



^logf^ 





1/2 






). 




/ ^oo 








Jb 





> log i /2 v^ 

■0oo V'oo 



1/2 



< |Vn"| 2 + 



(JV 



n\2 



^logf^ 

B \ Woo 



B 



and using the a priori bound (I42p . we see that the last term in uniformly bounded in L 1 ((0, T) x 

n). 

To bound the second term in (|73j) . we first use the inequality xy < C(x 2 log 2 (x) + ^ - ) 
for x,y >2 and then apply Jensen inequality. Hence, 
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(77) 
(78) 

(79) 



Vu 11 



log 



1p n Ipo 



< IWf + 



B ijjoo 1 - \R\ 

i r r v $ 



dR 



< Wu n \ 2 



+ 



1 



dR 



- 2 




log 2 





(N 



n\2 



4> n 



v' 



— log — 



V'oo 1 
2 



\R\ 



dR 



B roo 



V'c 



1 



l-RI 



and the last term can be bounded as in (|74[) . We notice here that the last inequality implies 
in particular that |r n,<5 ' K | 2 is equi-integrable in L l for fixed 6 and k. This is actually a very 
important fact that will be used again later. 

Proof of Lemma l5.2l To prove this lemma, we use dominated convergence and monotone 
convergence. Indeed, |2: ri ' <5 ' K | is decreasing in S,k, namely for < 5 < 5' and < k < k' , we 
have 



(80) \z n ' 5 '' K '\ < \z n ' S ' K \ < \Vu n W 

Passing to the limit weakly in n, we deduce that 



(81) \ z n,S'^'\ < \ z n,8,K\ < |Vu n (^ r 



-S,i 



and by monotone convergence, we deduce that G = \z n > s > K \ ' £ X and that for all < 5 and 
< k, we have \z S ' K \ < G. Moreover, we have 



(82) \z 5 ' K -z 5 '' K '\< 



y n,8,K\ \^n,5',K'\ 



Hence, there exists g G X such that z S ' K converges strongly to g in X. Now, we would 
like to prove that the limit g is equal to Vu n ip n . This follows from the fact that Vu n ijj n 
is equi-integrable in L 1 ((0, T) x f2 x B;dtdxdR) (without the weight). Indeed, denoting 
3>(x) = |x| log 1//2 (l + \x\), we have 

(83) $(|Vu> n |) < |WV n |(log(l + + log(l + \Vu n \)) 

(84) <^[|Vn"| 2 + log(^)] 

which is clearly bounded in ^(dtdxdR). 

5.3. Identification of J B 1 ,^3, dR. In this subsection, we give a relation between j3 and 

some defect measure related to the lack of strong convergence of Vu n in L 2 . To state the 
main proposition of this subsection, we introduce few notations. Let u n = v n + w n where v n 
and w n solve 



(85) 



(86) 



d t v n - Av n + VpJ = V.r n 
v n (t = 0) = 



d t w n - Aw n + Vp^ = -n n .Vu r ' 
u"(i = 0) = n n (t = 0). 



We further split w n into ivf + to £ where to™ is the solution with zero initial data and is 
the solution with zero right hand side. 
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In the rest of this subsection we will use d to denote 5, k. We define v n ' S = v n ' S ' K the 
solution of 

Q fV n,S _ Ay n,S + Vp n,6 _ y T n,S 
v n > S (t = 0) = 



(87) 



Extracting a subsequence, we assume that (T n ' S ,Vv n ' S , Vv n ,Vw n ) converges weakly in L 2 to 
some (t s , X7v s , Vv, Vw) and that 

(88) \X7v n > s \ 2 = |V/| 2 + / 

for some defect measure fi s £ M((0,T) x $7). We also denote /j, the limit of /j, s when 5 and 
then k go to zero (extracting a subsequence), namely 

(89) fi = lim lim = lim fi S . 
Proposition 5.4. We have 

(90) ^ = -j B ^l^jk\ dR - 

Proof of Proposition 15.41 

We introduce the following weak limits 

(91) T n > s : Vv n > s = W 55 

(92) T n ' s : Vv n = W 8 . 

Step 1: First, we would like to prove that W Sd and W s have the same limit W when d 
goes to zero and that this limit is in L , To prove this, we introduce for M > 0, the following 
weak limits 

(93) t™> 5 1| t „, 5 |< m : Vv n > s = Wff 

(94) T^ s l lrn , S]>M : Vv n > s = W ss - W 6 J 



M 



(95) T n ' S l [T u, Sl < M : Vv n = W t 
and 

(96) F^V^STF = Gm [^f = G s . 

Since for a fixed S, \r n ' S \ 2 is equi-integrable, we deduce that G S M converges to G 5 in L 1 when 
M goes to infinity and is monotone in M. Also, by monotone convergence, we deduce that 
there exists G € L 1 such that G s converges to G in L 1 when d goes to zero. Actually, G is 
the weak limit of |t"| 2 in the sense of Chacon. 

Let us fix e > 0. We choose 3$ and Mq such that for S < So and M > Mq, we have 
\\G - G s \\ L i + \\G - G s M \\ L i < e. We have 



(97) |t"> 5 | 2 = |rMl |rrM| < M |2 + |r-.n |T „ i5|>M |2 

(98) =Gf 1+ (G S -G S M ). 

Hence, we deduce that for S < Sq and M > Mq, we have for all n, \\T n,s l\ T n,s\ >M \\ 2 L2 < e 
and hence, by Cauchy-Schwarz we deduce that ||W** — Wjf/H^i < Cyfe and that \\W S — 
WmIIl 1 < Cy/e. Hence to prove that lim<5 W ss = liuis W s , it is enough to prove it for the M 
approximation, namely that 

(99) limW$ = ]imW&. 
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To prove (|99p . we first notice that r n,s — r n goes to zero in L p for p < 2 when 6 goes 
to zero uniformly in n. Then, by parabolic regularity of the Stokes system, we deduce that 
HV-i/M 5 _ Vv n ||ip((o,T)xf2) g° es to zero when 8 goes to zero uniformly in n for p < 2. Hence, 
(199]) holds. 

Step 2: In this second step, we will compare the local energy identity of the weak limit of 
(I87p with the weak limit of the local energy identity of (187p . 

On one hand, passing to the limit in (187]) and multiplying by v s , we deduce that 



(100) dj-^- - A^- + |Vw*| 2 + div(pfv*) = div(/.r 5 ) - Vu* : t s 



On the other hand, reversing the order, we get 



\ v 8\2 \ V S\2 



(101) d t ^- - A-^— !- + \Vv s \ 2 + H5 + dw(p s lV d ) = div(v s .t s ) - W ss 



For a justification of these two calculations, we refer to [38]. Comparing (|100p and (jlOip . we 
deduce that W 65 = Vv d : t s -fig- 

We would like now to send S to zero. 
First, it is clear that t s converges strongly to r in L 2 when d goes to zero. Hence, Vu* 5 also 
converges to X7v in L 2 . Besides, from the energy estimate, we recall that u n is bounded in 
L°°((0, T); L 2 (ft))nL 2 ((0, T); H 1 (U)) and hence by Sobolev embeddings that u n is bounded in 

2(D+2) D+2 

L d ((0, T) x SI) and that u n Vu n is bounded in ((0, T) x Si). By parabolic regularity 
of the Stokes operator applied to (|86p with zero initial data, we deduce that Vw" is bounded 

D + 2 -. D+2 D+2 

in Ld+i ((0,T); W 'wQ) and that dtw™ is bounded in L^+i ((0,T) x Si). Since r n is bounded 
in L 2 , we deduce from (|85p that Vf n is also bounded in L 2 ((0, T) x 0) and hence Vro" is also 
bounded in L 2 ((0,T) x SI). Moreover, it is clear that Vro^ is compact in L 2 ((0,T) x SI) and 
hence Vw™ is also bounded in I? and from the previous bounds on Vw™, we deduce that Vw" 

= — 8, K 

is compact in L p ((0,T) x SI) for p < 2. Hence, we deduce that Vw n : r(V> n ) = Vw : r(V ; ) 
(where we have used that r n,s is equi-integrable for each fixed 8) and from Corollary 15.31 that 

lima W ss = lim 5 W 5 = Vv n : r(^ n ) = Vv : r(^) + f B fcj-^vjfctdR. Finally, we deduce that 



H = lima^o V s = ~ J B N rrm dR - 



W\' 



5.4. Gronwall along the characteristics. Taking the difference between (j60|) and (|54p and 

dividing by iV 4 , we get (to be more precise, we have to take the difference between (|59p and 
(|54p . divide by A^ 4 and then send d to zero): 
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1 

W 



Vn":r(f) + Vtt : r ( Voo (V - 2^0'(^-)J 



(102) 



iV 4 
2ak 



4|Vi?/"| -2|V R <7| 2 + 

■5,K 



|V ii /«| 2 (log 1 /2 + iog-3/2 )( |!L ) 



2.9 



/^v«-(bg(£) + i) -(2e'(£)v^),) T 

y V>oc [7u^V 5 - V^Wii : (7 - l')g] 



RiRj 

Hip 



2 

~ iV 4 " 

4 

i=i 

where we denote the 4 terms appearing on the right hand side by Ai , 1 < i < 4 and we 
also denote 77 = iV™ ' K — f B g 2 ipoo = Jb[(9 h ) 2 ~ S^ipoodR. It measures the lack of strong 
convergence of g n to g in L 2 (dtdxtpoodR) . Notice that by the choice of the normalizing factor 
N, the defect measure is in L°°. 

First, we prove that A2 is nonnegative, namely we have the following lemma 



Lemma 5.5. We have 
(103) 

for some constant c. 

For the proof, we rewrite \Vng\ 2 as 



V ii r(log 1 / 2 (/n)2) + V ii r(log- 1 / 2 (/n)2) 



(104) |V R5 | 2 = V fl /" (log 1/2 (/™) 2 + log- 1 / 2 (/" )2) 

(105) 

(106) + 2V fl /» (log 1/2 (/ n ) 2 ) • V fl /» (log- 1 /2 (/ n )2) . 
Hence, we deduce that 

(107) A 2 = -^J^ 00 ( a + /3 + 1 ) 
where a, (3 and 7 are given by 



(108) 



a |V^I 2 0og 1/2 (£)) r , locl/2f ^. 
2 " /» J ^ Voo' 



(vr)io g i / 2 ( 



^00 



(109) 



g _ |V*/1 2 (log- 3/2 (fc)) 
2 /" 



/ nl og V2 ( ^ ) _ (V/ n )log -l/2 ( ^ ) 
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(110) 1 = 2 |vr|2 5 ' K -2(Vr)log 1 / 2 (^) (V/^log- 1 /^) 

2 Ipoo Woo 

We introduce the Young measure Vt t x,R,(A, A) associated to the sequence (Vf n , f n ) where 

n 2^' K 

A G M and A£f. Hence, the defect measure |Vi?(/ n — /)| satisfies : 

(111) |V i? (/--/)| 2 >|V R (r-/)| 25 ' K > /|A-| AV M , fl (A',A')| 2 ^, R (A,A) 

(112) =^ y'|A-A'| 2 ^ XiJi (A',A')^ >a!> ji(A,A) 

2<5> K 

Indeed, it is easy to see that |V#(/ n — /)| is bounded from above by the weak limit and 
from below by the Chacon limit of |V#(/ n — f)\ 2 - In the sequel, we will drop the t,x and R 
dependence of v and will denote v' = u(A', A') and v = u(A, A). Besides, a,j3 and 7 satisfy 

(113) a> J J A(A,\,A',\')v(A',\')v{A,X) 

and the same for (3 and 7 with A replaced by B or C where A, B and C are given by 

(114) 

A = l A | 2l0 g 1/2 ( A2 ) A / l0g l/2 ( y )2 + IA? W /2 (A0 2 AlQgl/2(A2) _ 2A . A / log l/2 (A 2 )log l/2 (A0 2 

A A 

(115) 

B = ^Ba'Io^(A') 2 + I A/ 1 2 log" 372 ( A/ ) 2 A 1 Qg i/2 ( A 2 j _ 2A.A'log- 1 /2 (A 2 )log -i/2 ( y )2 
A A 

(116) 

C = 2|A| 2 + 2|Af - 2A.A'( log 1 / 2 (A 2 ) log" 1 / 2 ^' 2 ) + log" 1 / 2 (A 2 ) log 1 / 2 (A / ) 2 ' 



To prove lemma [5T5j it is enough to show that A + B + C > §|A — A'| . First, we rewrite 
A + 5 + C as 

(117) A + B + C= \A\ 2 B 1 + (A'l 2 ^ - 2A.A'B 3 

(118) = |A - A'| 2 + lA) 2 ^! - 1) + |Af (S 2 - 1) - 2A.A , (5 3 - 1) 
where B\ , B2 and B3 are given by 



A Alog d/z (A 2 ) 



( m) B3 - ^ M ^ + logl/2(A2) ; ogl/2(A , 2) + + 

Actually, we will prove that if a is chosen big enough then (B\ — 1)(-B 2 — 1) > (-B3 — l) 2 
from which we deduce that A + B + C > |A — A'| 2 and the lemma would follow. Indeed, after 
simple calculations, we get 
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(B 1 -l)(B 2 -l)-(B 3 -lf 
Iog 1 / 2 (A 2 )Iog 1 / 2 (A' 2 



A V g log 1/2 (A 2 ) log 1 / 2 ^ 2 ) 

A' A log 1 /2( A '2) log 1 /2( A 2) 



+ 2 



+ 



logV 2 (A 2 ) + log 1 / 2 (A' 2 ) 



log 1 / 2 (A' 2 ) log 1 / 2 ^ 2 ) 



log 1 / 2 (A 2 )log 1 / 2 (A' 2 ) 



Alog(A 2 ) | A'log(A' 2 ) | 2 2 log 1/2 (A 2 ) 2 log 1 / 2 (A^ 2 ) 



A'log(A' 2 ) Alog(A 2 ) 



log 1 / 2 (A /2 ) log 1 / 2 (A 2 ) 



We will prove that the three terms appearing inside the brackets are nonnegative. This is 
obvious for the second one since it is of the form x + ^ — 2 for some x > 0. We recall that 
since (/ n ) 2 > a, we get that A > y/a on the support of v. For the first bracket, we assume 
that A' > A and write A' = A(l + e). Hence, the term in the first bracket is given by 



(122) 



1 + £+ 1 +2 _ 2 , 1 + MI±£i_ 2 . ' 



1 + E 



log A 



1 + 



log(l+e) 
log A 



and one can check easily that if A > \fa is big enough then (|122p is nonnegative. The same 
argument can be used for the third bracket. This end the proof of lemma [531 

To bound A%, we first observe that 
(123) 

and hence, 



2 - 2^G'(^) = - / nlog 1 / 2 (/n)2 /n lo g-l/2 (/ n)2 



(124) A 1 

(125) 
(126) 



1 

~w 

1 

"iV 4 

i 

w 



-<5,K 



Vu n . T (^n ' + v-u : r(V> - fa(/"log 1/2 (/ n ) 2 /"log~ 1/2 (/") 2 )) 



Vu : r(V - fa(/" log 1/2 (/") 2 /" log" 1 / 2 ^) 2 )) 



By convexity, it is clear that (/" - f) 2 = {f n f - f 2 > (/ n ) 2 - f n log 1/2 (/™) 2 /" log" 1/2 (/") 2 
and hence, 

(127) 

Voo(/" log^Cf 1 ) 2 /"fog- 1 / 2 !/") 2 ))! < (7 fa(/"-/) 2 / <Uv(/™-/)| 2 ^ ' 

vie ie 

Hence, 

- 

(128) 



(129) 



- -r I A'' 4 10iV 4 7 B 



fa|V(/"-/)| 2 
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The term between parentheses in the definition of A3 can be written as 



(130) j^- (logger) 2 +io g - 1 / 2 cr) 2 ) [/"log 1 / 2 !/™) 2 - /"iog 1/2 (/") 2 

If we denote ^^(II, A) the Young measure associated to the sequence (V x u n ,f n ), then 
we see easily that ^3 is given by 



( 



i3i) A z = -wJJ J (T (logl/2 A ' + log_1/2 a2) - ^ (logl/2 A ' 2 + log_1/2 A ' 2; 



(132) (Alog l/2 A 2_ A / log l/2 A /2 ) _R^^ oo dvdv f dK 

1 — \R\ Z 

The absolute value of the two factors inside the integral can be bounded respectively by 

log 1/2 A 2 log 1 / 2 A' 2 log 1 / 2 A 2 log 1 / 2 A' 2 

(133) in - n'|(^__^_ + ^_^ } + (|n| + |r[/|)( log__^ _ and 

(134) |A-A , |(log 1/2 A 2 +log 1/2 A /2 ). Hence 



(135) \A,\< — A I I hn-^^^ + ^^f—^—^dvdv'dR 
\ 10A^ 4 JbJ J A A' ; 1 - |i?| 2 y 

(136) + wJ B I /(l n l + l n 'l)l A - A, | 2 ^ + ^A#) 2 Y^^dudu'dR 

C f /"/"/..,„,. , . ,| 2 1 



(137) + N*j B j J (1 + |n| + |II/|)|A ~ V| 1 - \R\^ °° dudV ' dR 

1 1 C 

(138) < tuH jkH 7 \Vu\ 2 r] 

V ; ~ lOiV 4 ^ 10iV 4 iV 4 11 



Finally, to bound — A±, we split it into two terms : 
(139) \A\\ <-LJ ^ 00 \ lij \\Vg\ dR 



(140) < y^IVu"- Vn|| 5 ' K + ^^( 5 «- 5 )|V^|^ 

(141) j B ^ R 9\ 2 j^JF^W 

To bound A\, we first consider the case k > 1 where the term can be treated as A] usine 
1.2, ^ 2 f ... , g 



Ml 







>(l7ij 


< 






1 


C 


< 


10iV 4 ^ + 


w u 




1 


£. 


< 


10AT4^ + 


AT4 



dR 

\9? 



(142) < + — A I ^(gn - g y dR / ^00 |)2 dfl 



In the case fc < 1, we have to use (fTUj) instead of (J7J). Take < /3 < k and 7 = 
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(143) 



ipooS tIt^t dR 



2 f (g n -5) 2 log" 7 (/ n ) 2 

< / V>oo , 4 dR 



/ g 2l °g 7 (/") 2 Jn 

B ^°°(l-|i?|2)l+/3 °^ 



, B '(i-W) 7 ~^ ro ° (i - l^l 2 ) 1 "' 3 

To bound the second term, we use the following Young's inequality for a, b > 1 ab < 
i 

a log 7 a + e^ 7 -*. We denote d = 1 — |-R| 2 and hence 



g 2 log y (f n ) 2 /" 



n 2 r? 
9 j 7 ^ , r „,2 



di2. 



On the set {g 2 > ^} where e = we have log 7 -S+p < Clog 7 g 2 . Besides, we have using 



5 2 log 7 5 2 , D , / / >, 
Woo — rrzfl — dR < / V'ooS' log g 



1-/3 
2 



V'oo(|Vk 5 | 2 + (7 2 



B 



1+/3 
2 



'B ^ VB 

On the set {g 2 < ^}, we have 

which is integrable in the ball B with the measure ipoodR. 

To bound the first term on the right hand side of (j!43p . we first notice that 



(g n - g) 2 log- 7 (/«)2 < (fn - f)2 log 1 - 7 (C + (/" - f) 2 ) 
which can be easily proved using Young measures. Besides, we have using (|10p 



B 



< 



(i- \R\ 2 y-p 1 

^oo(r-/) 2 iog(c+(r-/) 2 ) 



B 



1+1 

2 



^o|VkCP-/)| ; 



B 



<-Srf/ ^oo(/ n -/) 2 iog(c + (r-/) 2 )) ( / ( ^|v«(/"-/)|- 



B 



for each A > 0. Passing to the limit weakly (more precisely, applying F n ' K ) to both sides and 
optimizing in A, we deduce that, 



(144) 



-g\g , D \ . c 

W) dR ) -w 



ipoog 2 log 



1-/3 
2 



N^Sb^-W-) J / VB 

Putting all these estimates together, we deduce that 



V>oc(|Vfl<?r + 5 



1+/3 

2 , 2n 1 2 i±£ i: 



7/ 2 ZU 2 



(145) 



; 2 i 2\ 1+yS 

V'ooff log 5 ry. 



B 



We can take /? — 0. Next, we observe that V'oo 

g 2 log g 2 dR < CN 2 . Indeed, if we 
introduce h n = g n log 1 / 2 g n , we see that N£ > (f B ipoo(h n ) 2 ) 1 ^ 2 and then it is easy to see 
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using that (x, y) — > 2- is convex that 



from which we deduce the claim. Hence ()145|) becomes 



1/2 



(146) 

±^(t,X(t,x)) + ^(t,X(t,x)) < C \Vu\ 2 ^L(t,X(t,x)) + C7 



1 + 



2V 



N 



First notice that the right hand side of (|146|) is in i 1 ((0, T) x K) for any bounded measurable 
set of (To prove this, we can observe that ^ is bounded and that using (|47|) . the term 
between brackets in (|146p is in L 1 ((0,T) x if)). Hence (fll6]l is well justified in the sense of 
distribution. In particular this justifies all the calculations done in this subsection starting 
from (fT02]l . 

Now, since the term between brackets in (|146|) is in L 1 ((0,T) x K), for almost all x, 



Jo l + J B V>c 



|v fl g| 2 

1 N 



(t,X(t,x)) is finite. Besides, for almost all x, N(t,X(t,x)) (which is 



constant in t) is bounded. Hence, we deduce that for almost all x, N 3 1 + J B ?/>oo 



+ 

| Vu| 2 (i, X(t, x)) is finite. Hence, by Gronwall lemma, we deduce that for a.e x, we have for 
all t < T, ^4 < epr(.x) an j s i nce ^(0, x) = due to the initial strong convergence, we 

deduce that v ^ x ^ = and hence r\ = and we deduce the strong convergence of g n to g. 
This yields that (u,ib) is a weak solution of ([T]) with the initial data (uq,^). 

6. Approximate system 

In the previous section, we proved the weak compactness of a sequence of solutions to the 
system (pQ). Of course one has to construct a sequence of (approximate) weak solutions to 
which we can apply the strategy of the previous sections. The only thing we have to make sure 
is that the calculations done in the previous section can be made on the approximate system. 
We consider a sequence of global smooth solutions (u n , ip n ) to the following regularized system 
where k is some integer that depends on D. In particular one can take k = 1 for D = 2 or 3: 

' d t u n + (u n ■ V)u n - vAu n + ±(A) 2fc u™ + Vp n = divr n , divu = 0, 



(147) 



d t ip n + u n .Vib n = div R 



^ T % = J B (R i ®V j U)ip n (t,x,R)dR (VUip n + f3Vip n ).n = on dB(0,R ). 
with a smooth initial condition (uq,^q) such that (v,q, i/Jq) converges strongly to (uQ,t/jo) in 
L 2 (U) xi^fixB) and log -tpo +Pq^oo converges strongly to Vo log -ipo+Po^oo 
in L 1 (QxS). We also assume that ([6]) holds uniformly in n. In the case Q is a bounded domain 
of M. D , we also impose the following boundary condition u n = Au n = ... = (A) 2fc_1 u n = at 
the boundary dfl. 

We do not detail the proof of existence for the system (|147p . We only mention that we have 
to combine classical results about strong solutions to Navier-Stokes system with the study of 
the linear Fokker-Planck equation (see [51]). In particular the following operator was used 



(148) 



Lib = -div(ipooV - — ) 

Woe 
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on the space % = L 2 (-^-) and with domain 

(149) D(L) = Uenl^V^ en, divtyooV^en andVooV-^| 9B 

I *P°o Woo Woo 

Also the following two Hilbert spaces Ti 1 and Ti 2 are used in the construction : 

(150) H 1 = { ip G U | fipc 



■000 



2 

+ — dR < oo 



2 



,15!) W 2 = |v € W 1 /(**vi)) £<».} 

Following the the proof of existence given in |51| . we can prove 

Proposition 6.1. Take u% G ff s (fi) and > suc/i ^ Vo ~ Po^oo G # s (ty L 2 {^)) with 
Pq = J ifrydR G L°°(0). Then, t/iere exists a global unique solution (u n ,ip n ) to \147\ ) such that 
{u n ,^ n - p n Voo) «s in C([0,T);# s ) x C([0, T); fPQR*; L 2 {^))) for all < T. Moreover, 
u n G L 2 ([0,T);F s+A; ) and ^ n - p n Voo G L 2 ([0,T); ff'flR^;?* 1 )). 

Remark 6.2. T/ie proof is exactly the same as the proof of Theorem 2.1 of [51] with few 
differences: 

• In [51], we on/y /iad Zoca/ existence whereas here, we have global existence since we 
have more regularity. 

• Theorem 2.1 of [51 1 was stated in the whole space. Of course in the case of a bounded 
domain, we have to use energy bounds for Navier-Stokes written in a bounded domain. 

• In theorem 2.1 o/[5T] we assumed that J ip^dR = 1. The result can be easily extended to 
this more general case. We also point out that there is a small mistake in the statement 
of the theorem 2.1 of [51J. Indeed, one has to read if)Q — ipoo G H S ({1; L 2 (^-)) instead 

q/-0o G H s {yi;L 2 {-^p^)) when the problem is in the whole space. 

It is clear that the solutions constructed in Proposition 16.11 satisfy the free-energy bound 
(|32p and the extra bound (142p (with f2 replaced by K in the whole space case). 

Once we have our sequence of regular approximate solutions, we have to check that all the 
computations performed in the previous section can be done on this sequence {u n ,ip n ). The 
only point to be checked is that Proposition 15.41 still holds since the rest of the proof only 
involves the transport equation. Now, v n and w n solve 

d t v n - Av n + ±A 2k v n + Vpi = V.r n 
v n (t = 0) = n 



(152) 
(153) 



(154) 



d t w n - Aw n + ±A 2k w n + Vp% = -u n .Vu n 
v n (t = 0) = u n (t = 0) 

and we define v n ' S the solution of 

Q tV n,S _ Av n,S + l A 2k y n,S + Vp n,S = y T n,S 
V n > S (t = 0) = 

Step 1 of the proof of Proposition 15.41 is the same with the difference that one has to apply 
parabolic regularity for the perturbed Stokes operator which yields the same uniform in n 
estimate. Hence, we deduce that ||Vv n ' 5 — Vu n ||j>((o,T)xn) g° es to zero when S goes to zero 
uniformly in n for p < 2. 

For the second step, we first notice that (llOOp remains the same since ^{A) 2k u n converges 
weakly to zero. Moreover, multiplying the first equation of (|152p by v n,s , we get 
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L.n,<5|2 \„.n,S\2 

(155) dt^-^- - + \Vv n ' S \ 2 + G n ' 5 + divO^'V'*) = div(v n > s .r n ' d ) - r"' 6 : Vv n ' S . 

where G n,s is given by 



G n,S = I[di Vi (ViA 2fe - 1 ?; n ' <5 .?; n ' <5 - A 2k ^v n ' d .V^™' <5 + 

A k v n > s .V i A k - 1 v n < 5 ) + A k v n ' S .A k v n > 6 } 



(l^v n A 2t-2.,ni a„,«,<5 A L.n J n A fe— l„,ra 



Using the fact that ± f Q \A k v n > s \ 2 and J c , / n |Vv n > 5 | 2 are uniformly bounded, we deduce 
easily that G n ' S = \A k v n ' S \ 2 > and hence passing to the limit in (I157p . we deduce that 

(157) d t ^Y~ - + |V/| 2 + n s + dw(pfv s ) < div(v s .T S ) - W ss . 

and hence Proposition 15.41 is replaced by an inequality fi < — J B fyj | dR which is the 
inequality that we need in the rest of the proof. 

7. Conclusion 

In this paper we gave a proof of existence of weak solutions to the system (TjQ), using the 
fact that a sequence of regular solutions to the approximate system (|147j) converges weakly 
to a weak solution of (pQ). We would like here to mention few important open problems (with 
increasing level of difficulty, at least this is what the author thinks): 

• The zero diffusion limit in x. If we add a diffusion term -A. x if} in the Fokker-Planck 
equation of (H|), then one can prove the global existence of weak solutions to the 
regularized model. A natural question is whether we recover a weak solution of the 
unregularized system (pQ) when n goes to zero. This is the object of a forthcoming 
paper [53] • The difficulty comes from the fact that the calculation of section [5] used 
in a critical way the fact that we had a transport equation in the x variable. 

• Relaxing the assumption ()42|) . This extra bound was only used to give some extra 
control on the stress tensor. Can we prove the same existence result without it ? 

• Other models. A natural question is whether we can extend this to the Hooke model 
(where the system can be reduced to a macroscopic model). We were not able to 
perform this. The main difficulty is that we do not know whether the extra stress 
tensor r is in L 2 . Nevertheless, we know how to use the strategy to this paper to 
prove global existence for the FENE-P model |52j . 

• Regularity in 2D. Many works on polymeric flows are motivated by similar known 
results for the Navier-Stokes system. In particular a natural question is whether one 
can prove global existence of smooth solutions to (pQ) in 2D. We point out that this is 
known for the co-rotational model \47\ [51] . Of course this seems to be a very difficult 
problem since, we only have an L 2 bound on r and that an L°° bound on r was 
necessary in the previously mentioned works. In particular the similar result is not 
known for the co-rotational Oldroyd-B model where one can prove IP bounds on r for 
each p > 1. 

• Is system (pQ) better behaved than Navier-Stokes. One does not expect to prove results 
on (pQ) which are not known for Navier-Stokes since (pQ) is more complicated than 
Navier-Stokes. However, one can speculate that due to the polymers and the extra 
stress tensor, system (pQ) may behave better than Navier-Stokes and that one can 
prove global existence of smooth solutions to (pQ) even if such result is not proved or 
disproved for the Navier-Stokes system. 
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